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James  M  Rondinelli 

Department  of  Materials  Science  &  Engineering,  Northwestern  University,  Evanston,  Illinois  60208,  USA 


I.  Program  Summary 

The  program  utilized  electronic  structure  based  methods  to  de¬ 
sign  noncentrosymmetric  metals  (NCSM)  in  transitional  metal 
oxides  (Figure  1).  The  report  describes  the  progress  made  to¬ 
wards  understanding  why  any  NCS  metals  should  exist,  how  to 
find  them,  and  the  physical  properties  of  these  materials  with 
potential  applications. 

Main  outcomes  of  the  program  include  (a)  establishing  a 
weak  electron-lattice  coupling  principle  as  a  platform  for  ex¬ 
plaining  the  stability  of  NCSM  phases,  (b)  formulating  a  tax¬ 
onomy  and  scheme  to  aid  in  the  design  of  new  NCSM,  (c)  de¬ 
signing/predicting  new  metallic  transition  metal  oxides  without 
inversion  symmetry,  and  (d)  efforts  by  experimental  colleagues 
in  realizing  the  NCSM  predictions. 

The  key  results  include  the  identification  of  new  layered 
ruthenates,  cuprates,  osmates,  and  molybdates  that  fulfill  the 
formulated  operational  principle.  The  culmination  of  this  effort 


Figure  1:  Metals  with  noncentrosymmetric  crystal 
structures  are  rare,  but  this  class  of  compounds 
with  finite  band  occupancy  and  broken  inversion 
symmetry  (inset)  may  have  desirable  physical 
properties  for  device  applications.  This  program 
delivered  a  design  framework  for  predicting  new 
noncentrosymmetric  conductors  and  proposed 
new  technology  arenas  where  they  may  be  utilized. 

has  also  spawned  numerous  acentrix  oxides  and  ferroic  discoveries  from  mechanisms  that  emerge  in  layered 
oxide  compounds  otherwise  not  possible  in  simpler  structures  and  chemistries. 

_  (a)  _  _  _  (b)  _  _  _(c) 

II.  Scientific  Discoveries  &  Accomplishments 

A.  New  Oxides  without  Inversion  Symmetry. 

Representation  theory  analyses  performed  during  this 
program  examined  how  changes  in  the  direction  of 
cation  ordering  and  octahedral  rotations  can  lead  to 
the  loss  of  inversion  symmetry  in  A I  A'  III  superlattices: 
layered  (OOl)-type,  columnar  (llO)-type,  and  rock-salt 
(lll)-type  ordering  (Figure  2).  In  this  case,  the  superposition  of  the  cation  ordering  and  rotations  lead  to 
31  unique  space  groups  and  11  unique  crystal  classes.  Six  of  these  crystal  classes  are  noncentrosymmetric. 
The  important  result  is  that  there  are  no  combinations  of  octahedral  rotations  and  columnar  ordering  that 
produce  noncentrosymmetric  (NCS)  crystals,  sans  any  other  distortions.  We  identified  that  the  most  promising 
compositions  to  search  for  acentricity  in  A-site  layered  or  rock-salt  ordered  1/1  superlattices;  we  reported  these 
results  in  an  article  published  in  Chemistry  of  Materials.  This  work  was  also  extended  to  include  B/Br,  and 
mixed  A/A',B/B'  1/1  superlattices  (Figure  3).  This  structural  ‘heat  map’  establishes  a  design  tool  that  makes  it 
possible  to  select  chemistries  and  polymorphs  in  perovskite  oxides  that  are  conducive  to  NCSM  behavior. 

Our  results  reveal  that  the  ordered  arrangement  of  cations  plays  a  critical  role  in  determining  the  allowed 
rotation  patterns  that  can  lift  inversion  symmetry.  The  main  conclusion  is  that  rotations  of  oxygen  octahedra 


Figure  2:  The  three  types  of  A/A'  cation  ordering  in¬ 
vestigated:  (a)  layered,  (b),  columnar,  and  (c)  rock  salt 
arrangements.  Oxygen  atoms  are  omitted  for  clarity. 
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Figure  3:  The  type  of  space  group  symmetry  produced  by  combining  each  of  the  layered  ([001]),  rock  salt  ([1 1 1]),  and  columnar 
([1 1 0])  type  A-  and  B-site  cation  ordering  patterns  with  the  1 5  Glazer  rotational  patterns  of  the  B06  octahedra.  540  distinct 
structures  are  represented,  produced  by  also  considering  [100]  and  [010]  layered  ordering,  and  [101]  and  [01 1]  columnar  ordering. 
Each  separated  block  depicts  one  ordering-rotation  combination;  these  are  further  subdivided  into  whether  a  centrosymmetric, 
polar,  enantiomorphic  (chiral),  or  SHG-active  space  group  is  exhibited  (indicated  by  a  filled  square)  by  the  given  combination. 
The  number  of  distinct  space  groups  produced  is  given  by  the  number  in  each  of  the  filled  squares.  Note  that  because  many 
non-centrosymmetric  space  groups  fall  into  more  than  one  category  (for  example,  point  group  2  is  both  chiral  and  polar),  the 
numbers  do  not  sum  to  the  total  number  of  structures. 


can  induce  ferroelectric  polarizations  if  the  cations  are  arranged  in  a  layered  or  rock  salted  order.  Using  this 
chemistry-agnostic  understanding,  we  then  designed  several  new  polar  oxides.  Our  main  discoveries  include: 


o  New  Chiral  Perovskites.  We  studied  rhombohedral  instead  of  orthorhombic  perovskites  in  superlattices 
of  rare-earth  aluminates.  We  found  that  although  most  of  the  compounds  are  polar  and  exhibit  a 
small  spontaneous  polarization  of  2  uC/cnB,  these  compounds  undergo  a  transition  to  a  chiral  phase 
accompanied  by  a  dramatic  increase  in  the  piezoelectric  response  that  persists  to  high  temperature 
(Figure  4),  which  could  inspire  new  high  temperature,  low-fatigue  piezoelectric  devices.”  The  main 
finding  is  that  while  only  (La,Pr)Al20e  and  (CeyPfiAfiOe  are  polar,  all  were  found  to  undergo  a  transition 
to  an  R32  chiral  phase  at  higher  temperatures.  By  computing  the  total  polarization,  we  found  that  chiral 
(La,Nd)Al206  has  zero  polarization  while  the  polar  phases  have  values  of  ~2pC/cm2.  This  rare  example 
of  chiral  perovskites  was  reported  in  Physical  Review  B. 


o  New  Ferroic  Transition — Per  r /electricity.  We  identify  a  first-order,  isosymmetric  transition  between 
a  ferrielectric  (FiE)  and  ferroelectric  (FE)  state  in  A-site  ordered  LaScCVBiScOs  and  LalnOs/BilnOs 
superlattices.  Such  a  previously  unreported  ferroic  transition  is  driven  by  the  easy  switching  of  cation 
displacements  without  changing  the  overall  polarization  direction  or  crystallographic  symmetry.  Epitaxial 
strains  less  than  2%  were  predicted  to  be  sufficient  to  traverse  the  phase  boundary,  across  which  we 
captured  a  ~  5X  increase  in  electric  polarization.  Unlike  conventional  Pb-based  perovskite  ceramics 
with  a  morphotropic  phase  boundary  (MPB)  that  show  polarization  rotation,  we  predicted  an  electrome¬ 
chanical  response  up  to  102  pC/N  in  the  vicinity  of  FiE-FE  phase  boundary  due  to  polarization  switching 
without  the  change  in  symmetry.  The  significance  of  this  result  is  that  the  transition  is  an  alternative 
ferroic  transition  to  obtain  sizeable  piezoelectric  responses,  with  the  additional  advantage  of  occurring  at 
integral  stoichiometry  in  an  ordered  superlattice,  without  introducing  compositional  complexity,  toxicity,  or 
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a  change  in  symmetry— features  absent  in  conventional  Pb-based  piezoelectrics.  This  work  was  published 
in  Advanced  Materials  Interfaces. 3 


o  Magnetic  Oxides  with  Enhanced  Polarizations.  Recent  experimental  and  theoretical  work  has  shown 
that  the  double  perovskite  NaLaMnWOf,  exhibits  antiferromagnetic  ordering  owing  to  the  Mn  d  states, 
and  computational  studies  further  predict  it  to  exhibit  a  spontaneous  electric  polarization  due  to  an 
improper  mechanism  for  ferroelectricity  [King  et  ai,  Phys.  Rev.  B,  2009,  79,  224428;  Fukushima  et  al, 
Phys.  Chem.  Chem.  Phys.,  2011, 13, 12186],  which  make  it  a  candidate  multiferroic  material.  Using  first- 
principles  density  functional  calculations,  we  investigated  nine  isostructural  and  isovalent  AA'MnWOf, 
double  perovskites  (A=Na,  K,  and  Rb;  A'- La,  Nd,  and  Y)  and  formulated  crystal-chemistry  guidelines 
describing  how  to  enhance  the  magnitude  of  the  electric  polarization  through  chemical  substitution  of 
the  A-site  while  retaining  long-range  magnetic  order.  The  simultaneous  presence  of  anti-ferromagnetic 
ordering  of  the  Mn  atoms  and  spontaneous  electric  polarization  makes  these  compounds  multiferroic, 
a  highly  desired  type  of  material  for  applications  such  as  four-state  memory  devices.  By  substituting  a 
variety  of  alkali  metal  (Na,  K,  Rb)  and  rare  earth  (La,  Nd,  Y)  elements  on  the  A-  and  A'-sites,  we  showed 
that  the  magnitude  of  the  ferroelectric  polarization  can  be  increased  by  up  to  150%  without  destroying  the 
magnetic  properties.  This  work  was  published  in  Dalton  Transactions. 


The  insulating  acentric  oxides  and  ferroic  mechanisms  described 
above  were  direct  outcomes  of  our  efforts  in  this  program  focused  on 
the  design  and  understanding  of  noncentrosymmetric  metals  (NCSM), 
which  are  described  next. 

B.  Noncentrosymmetric  Metals. 

B.l.  New  Classification  Scheme  for  Materials  Discovery. 

We  established  a  framework  to  design/discover  NCSM  oxides, 

which  originates  from  Anderson’s  work  on  “ferroelectric  metals,” 
where  he  writes  that,  “while  free  electrons  screen  out  the  electric  field 
[in  materials]  completely,  they  do  not  interact  very  strongly  with  the 
transverse  optical  phonons  and  the  Lorentz  local  fields  [that]  lead  to 
ferroelectricity,  since  umklapp  processes  are  forbidden  as  k  -*  0.”  We 
recasted  this  observation  into  an  operational  design  principle  that  states: 


Low  T  High  T 


Figure  4:  Predicted  piezoelectric  coeffi¬ 
cients  for  La^Pr^AICL. 


The  existence  of  any  NCSM  relies  on  weak  coupling  between  the  electrons  at  the  Fermi  level,  and  the 
(soft)  infrared-active  phonons  responsible  for  removing  inversion  symmetry. 


An  essential,  implicit,  materials  constraint  is  that  the  low-energy  electronic  structure  derives  from  an  electron 
count  giving  partial  band  occupation,  which  may  be  obtained  by  judicious  selection  of  the  cation  chemistries. 

After  surveying  families  of  known  NCSMs,  we  identified  three  routes  by  which  the  weak-coupling  op¬ 
eration  principle  may  be  validated  in  a  solid-state  system  (Table  1),  and  reported  our  findings  in  Nature 
Communications.  The  summary  of  our  comprehensive  survey  will  be  reported  later  in  2015/2016  in  the  form  of 
an  invited  article  to  J.  Phys.:  Condens.  Matter.  The  classification  scheme  entails  specify  the  symmetry  behavior 
of  the  lattice  phonons  that  would  characterize  the  inversion  symmetry  loss  across  a  centrosymmetric  to  NCS 
structural  transition  (Table  1).  Because  the  first  approach  is  ‘trivial’,  we  have  rather  focused  our  efforts  on 
exploring  NCSM  realized  by  Mechanism  2  and  3  in  the  layered  cuprates  (LnSr 2 C u 2 G a O 7 ,  Mechanism  2)  and 
ruthenate,  osmate,  and  molybdate  superlattices  (Mechanism  3). 


B.2.  Designed  Polar  Metals — Ruthenates. 

Using  our  operational  principle,  we  designed  ultrashort  period  SrRuOj/CaRuO;  superlattices.  Our  main 
results  reported  in  Nature  Communications 6  include: 
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Table  1:  Routes  to  lift  inversion  symmetry  in  metallic  compounds.  The  scheme  is  formulated  by  classifying  the  symmetry 
requirements  of  the  lattice  mode  instabilities  of  a  centrosymmetric  (reference)  phase  which  would  drive  a  thermodynamic 
transition  to  a  noncentrosymmetric  crystal  structure.  Although  compositional  ordering  may  be  an  obvious  route  to  lift  inversion 
symmetry,  we  recognize  that  experimentally  in  practice  it  may  be  challenging  to  achieve  it,  especially  in  oxide  materials  owing 
to  bond-coordination  requirements;  therefore,  although  the  route  is  simple,  it  is  by  no  means  straightforward.  Note  that  all 
materials  are  expected  to  show  some  interesting  physical  properties  enabled  by  an  acentric  crystal  structure  with  itinerant 
electrons,  independent  of  the  mechanism  leading  to  the  mutual  coexistence  of  the  prerequisite  atomic  and  electronic  structure. 

Inversion-lifting  Mechanism  Phonon  Requirements  Description 


(1)  Compositional  Order 
(*) 

(2)  Packing  of  Acentric 
Polyhedra 


(3)  Geometric-induced 
Displacements 


None  Achieved  by  decorating  one  or  more  interleaved  lattices  with  mul¬ 

tiple  cations,  e.g.,  as  realized  in  half-Heulser  alloys  or  tri-color, 
•/ABC/ ABC/-,  superlattices. 

1,  k  =  0  or  k  +  0  mode  Realized  by  the  alignment  of  acentric  metal-oxygen  polyhedral 
units  arising  from  either  out-of-center  metal  distortions  or  because 
they  intrinsically  lack  inversion  (MO4  tetrahedron),  and  contribute 
few  (if  any)  states  to  or  near  the  Fermi  level. 

>  2,  k  ±  0  coupled  modes  Obtained  by  anharmonic  coupling  of  two  or  more  centric  lattice 
modes,  which  cooperatively  remove  inversion  and  produce  cation 
displacements  that  do  not  gap  the  electronic  structure.  Accessible 
in  some,  effectively,  two-dimensional  compounds,  e.g.,  naturally 
layered  structures  or  3D  systems  with  bi-color  ordering. 


1.  The  ground  state  structure  is  polar,  space  group  Pmc2\  (Figure  5).  The  inversion  symmetry  is  due  to 
centric  RuOe  octahedral  rotations  that  induce  polar  Sr,  Ca,  and  remarkably  Ru  displacements. 


2.  The  electronic  structure  is  robustly  metallic  in  the  polar  ground  state.  Our  calculations  show  that  spin-orbit 
interactions,  which  play  a  significant  role  in  determining  the  orbital  structure  in  some  transition-metal 
oxides,  ’  weakly  modified  the  Ru  4 d- band  occupation  at  the  Fermi  level. 


3.  The  thermopower  coefficients  are  highly  anisotropic.  The  anisotropy,  which  we  quantify  by  computing  the 
contributions  along  the  principal  axes  and  then  taking  the  sum  ASL  =  Sc  -  ^  (Sa  +  Sb),  is  approximately 
6.3/iV/K  at  300  K  derived  from  its  polar  structure.  It  exceeds  that  of  YBa2Cu307_5  (~  7  -  10/iV/K).1" 


Our  study  on  the  SrRu03/CaRu03  1/1  su¬ 
perlattice  provides,  to  the  best  of  my  knowledge, 
the  first  examination  and  theoretical  calcu¬ 
lation  of  the  thermoelectric  properties  of  a 
NCSM.  The  implications  of  this  behavior  in¬ 
clude  the  potential  for  highly  anisotropic  thermo¬ 
electric  devices.  For  example,  the  peculiar  ther¬ 
mopower  anisotropy  in  the  SrRu03/CaRu03 
1/1  superlattice  dictates  that  the  electric  field 
resulting  from  an  applied  heat  flux  to  the  ma¬ 
terial  will  be  non-collinear.  This  property  is  a 
fundamental  feature  for  any  anisotropic  ther¬ 
moelectric  devices.  It  enables  the  heat  flux  to 
be  measured  in  a  geometry  perpendicular  to 


Figure  5:  Crystal  structure  of  the  layered  (Sr,Ca)Ru206  perovskite 
with  Ru06  octahedra  exhibiting  the  a+b~b~  octahedral  tilt  pattern  (a). 
The  out-of-center  Ru4+distortions  (bond  lengths  in  A)  produce  local 
dipoles  which  cooperatively  order  along  the  c-direction  (b),  removing 
the  mirror  plane  perpendicular  to  the  A-site  cation  order.  The  local 
coordination  environment  for  (c)  Sr  (tetra-capped  trigonal  prism)  and 
(d)  Ca  (bi-capped  trigonal  prism)  reveal  that  large  A-site  distortions  from 
the  ideal  12-fold  coordination  occur  in  opposite  directions  (arrows), 


along  the  polar  c-axis. 

the  induced  electrical  current,  specifically  at  locations  where  the  temperatures  are  equal  (Figure  6).  Given 
the  metallic  conductivity  of  the  ruthenate  superlattice  and  that  the  relaxation  time  is  small  in  comparison 
with  that  of  semiconductors  or  insulators  already  finding  use  in  thermoelectric  devices,  this  polar-NCSM 
would  enable  sensing  on  sub-nanosecond  timescales  (comparable  to  YBCO  ).  New  applications  of  these 
oxide  materials  could  be  found  in  ultrafast-thermoelectric  devices,  '  where  stability  under  extreme  conditions, 
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speed,  and  the  ability  to  measure  heat  fluxes  of  high  density  are  key  requirements,  e.g.,  thermal  (heat)  radiation 
detectors.  Moreover,  because  of  the  compatibility  of  perovskite  oxides  with  Si-based  CMOS  technologies, 
we  anticipate  these  designed  NCS  metals  will  more  readily  find  integration  and  device  development  than 
previously  identified  materials  with  large  anisotropic  thermoelectric  responses. 

Furthermore,  we  investigated  the  linear  op¬ 
tical  properties  of  this  superlattice  using  den¬ 
sity  functional  theory.  The  results  are  published 
in  /.  Phys.  Condens.  Matter.  26,  265501  (2014). 

In  this  effort,  we  found  that  the  polar  ferro¬ 
magnetic  metallic  ruthenate  exhibits  optical 
anisotropy  along  the  principal  directions  of  the 
optical  indicatrix  owing  to  the  absence  of  inver¬ 
sion  symmetry  in  the  crystal  structure.  We 
then  used  the  generalized  Drude  model 
and  deduced  an  enhancement  of  the  effective 
mass,  m*  ~  4.9me.  Upon  removal  of  the  polar 
distortions,  we  determined  there  was  a  decrease 
in  the  effective  mass  to  m*  ~4.4me. 

Our  main  discovery  is  uncovering  a  dependence  of  the  effective  mass  on  the  amplitude  of  the  polar 
distortion.  This  suggests  that  control  over  the  degree  of  inversion  symmetry  breaking  via  the  polar  displacements 
could  be  used  and  an  effective  ‘ handle ’  over  the  strength  of  the  electronic  correlations  in  the  ruthenate  superlattice, 
and  perhaps,  in  polar  metals  generally.  We  conjecture  that  epitaxial  strain,  which  has  been  successfully  used  to 
induce  (or  enhance)  polar  distortions  in  conventional  dielectrics  (ferroelectrics),  -  could  also  be  applied 
to  thin  films  of  SrRuO  j/CaRuO  ;  to  modulate  the  amplitude  of  the  polar  Sr  and  Ca  displacements;  thereby, 
tuning  the  electron-electron  interactions  which  would  lead  to  control  over  the  effective  carrier  masses. 

Lastly,  in  this  effort  we  have  and  con¬ 
tinue  to  work  with  Prof.  Chang-Beom 
Eom’s  Group  to  experimentally  realize 
the  3/3  period  ruthenate  superlattices 
(Figure  7),  which  is  more  amenable  to  ox¬ 
ide  thin  film  growth.  We  have  performed 
DFT  calculations  to  predict  the  ground 
state  structure,  the  electronic,  and  mag¬ 
netic  properties  of  this  superlattice,  and 
summarize  our  main  results  here:  we  im¬ 
posed  as  a  constraint  b  =  c  =  3.867  A,  sim¬ 
ulating  the  situation  of  thin -film  growth 
of  the  cuprate  on  NdGaOj.  The  ground 
state  structure  has  space  group  Pmc2\,  sim¬ 
ilar  to  the  1/1  superlattice,  and  exhibits  the 
a+b~b~  octahedral  tilt  pattern.  The  density 
of  states  is  consistent  with  that  reported 
in  Ref.  20  where  the  valence  band  is  com¬ 
posed  largely  of  O  2 p  states  hybridized  with  Ru  4 d  states,  with  the  oxygen  states  predominately  found  in 
regions  below  the  Fermi  energy.  The  Ru  atoms  are  magnetic  with  a  moment  of  0.154 pg  on  Ru(l),  0.4  ()()«/{ 
on  Ru(2),  and  0. 015  w  ^ .  These  superlattices  are  now  being  studied  using  second-harmonic  generation  by  the 
Gopalan  group  at  PennState  to  confirm  our  theoretical  predictions.  Results  from  this  work  is  planned  for 
publication  in  2015  or  2016. 


(Sr0.51Ca0  5)RuO3  /  NdGa03  (001) 


FWHM  ~  0.06° 


[(SrRuO3)3/(CaRuO3)3]10  /  NdGa03  (001) 


FWHM  -  0.06° 


20  O 


Figure  7:  Epitaxial  growth  of  a  (Sr,Ca)Ru03  film  and  a  [(SrRuO3)3/(CaRuO3)3]i0 
superlattice  on  orthorhombic  NdGa03  (001)  substrates  synthesized  in  the 
Eom's  Group  (University  of  Wisconsin,  Madison).  AFM  images  of  a 
(Sr,Ca)Ru03  film  (a)  and  a  [(SrRuO3)3/(CaRuO3)3]i0  (b)  superlattice  on  NdGa03 
(001)  substrates.  The  insets  represent  the  corresponding  RHEED  images 
of  the  CaRu03  and  SrRu03  films  (c),  6-26  scans  of  a  (Sr,Ca)Ru03  film  and  a 
[(SrRu03)3/(CaRu03)3]io  superlattice  on  NdGa03  (001)  substrates. 


Figure  6:  Thermoelectric  device  based  on  an  anisotropic  thermal 
conductor.  The  temperature  gradient  along  the  x2  direction  results  in 
an  electric  field  along  the  x\  direction,  which  can  generate  a  current 
density.  Note  that  along  the  x2  direction  there  is  no  temperature 
gradient,  i.e.,  dT /dxi  =  0,  and  thus  Ta  =  TB. 
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B.3.  Understanding  Conducting  Copper  Oxides  without  Inversion  Symmetry. 

We  studied  the  lattice  dynamics  and  electronic  properties  of  the  layered  copper  oxide  superconductor 
LaSr2Cu2Ga07  (LSCGO)  ’  as  a  model  system  to  explain  the  existence  of  NCSMs.  Our  main  conclusion  from 
these  investigations  is  that  the  symmetry  breaking  in  LSCGO  originates  from  the  chemistry  of  3d10  Ga, 
which  prefers  to  form  tetrahedral  coordinations.  The  incompatibility  between  displacements  that  derive  from 
the  zone-center  phonons  and  the  metallic  behavior  is  circumvented  in  LSCGO  by  the  fact  that  these  distortions 
are  largely  decoupled  from  the  electronic  structure  at  the  Fermi  level  (Figure  8) — they  occur  at  much  lower 
energy  and  hence  fulfill  our  weak-electron  lattice  coupling  principle.  We  showed  that  the  distortion  of  the 
Ga04  polyhedra  unit  leads  the  loss  of  inversion  symmetry  through  Mechanism  2. 

Owing  to  the  near  degeneracy  of  the  polar  and 
non-polar  LSCGO  low-symmetry  structures,  we  per¬ 
formed  a  series  of  calculations  to  understand  the  ori¬ 
gin  for  the  degeneracy  as  means  to  understand  how 
to  favor  the  polar  phase  over  the  non-polar  structure. 

Figure  9  depicts  the  Ima2  and  Pnma  structures.  The 
main  structure  difference  between  the  two  phases  is 
the  alignment  of  the  tetrahedral  Ga04  units.  In  the 
polar  Ima2  structure,  the  base-to-vertex  direction  of 
each  tetrahedra  is  aligned  along  the  same  (c)  direction 
[Figure  9,  left].  In  contrast,  the  non-polar  Pnma  struc¬ 
ture  possesses  GaC>4  that  are  anti-aligned.  The  different 
layers  point  either  along  the  +c  or  -c-direction. 

Remarkably,  the  electronic  structure  also  does  not 
depend  on  the  alignment  of  the  GaC>4  units  [Figure  9, 
right],  suggesting  to  us  that  the  LSCGO  electronic 
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Figure  8:  Electronic  properties  of  LaSoCibGaOy  for  the  centric 
Imma  (orange)  and  acentric  Ima2  structure.  Comparison 
between  the  total  DOS  (a),  Cu02  layers  (b),  Ga  d- states  (c),  Ga 
s-  and  p-states  (d). 


structure  is  more  accurately  described  by  two  “disconnected”  sub-blocks.  Figure  9  depicts  the  two  sub¬ 
blocks  (outlined  in  a  broken  line).  Each  block  contains  one  plane  of  Ga04  tetrahedra.  We  proposed  that 
because  the  inter-block  distance  d  is  large  under  the  equilibrium  constraints  id  =  3.57  A),  the  electronic 
structure  and  phase  stability  is  essential  independent  of  the  Ga04  orientation  (aligned  or  anti- aligned).  We 
tested  this  hypothesis  by  keeping  the  atomic  positions  of  the  structural  blocks  fixed  and  reduced  the  inter-block 
distance.  For  each  fixed  d- value,  we  then  computed  the  density-of-states.  We  find  that  the  blocks  become 


Figure  9:  Crystal  structure  for  the  polar  Imal  and  non-polar  (Pnma)  LSCGO  phases  (left).  The  main  difference  is  the  alignment 
of  the  Ga04  tetrahedral  units  (red  or  green)  along  the  crystallographic  c-axis.  The  molecular  sub-block  units  are  outlined.  Effect 
of  inter-block  spacing  d  on  the  electronic  structure  degeneracy  between  the  two  phases  (right).  Bringing  the  molecular  blocks 
closer  by  reducing  d  leads  to  interactions  that  split  the  electronic  equivalency  of  the  two  phases. 


j.m.rondinelli  |  jrondinelli@northwestern.edu 


6 


FINAL  REPORT 


(a) 


(b) 


“coupled”  when  d  is  approximately  2.16  A  as  seen  by  the  different  electronic  structure  between  the  lma2  and 
Pnma  phases  (Figure  9).  Based  on  these  experiments,  we  arrived  at  our  main  result:  The  polar  cuprate 
phase  can  be  favored  by  reducing  the  inter-block  spacing.  Our  subsequent  efforts  in  this  project  focused 
on  identifying  feasible  routes  to  control  that  space  and  thus  the  phase  stability. 

First,  we  examined  how  chemical  pressure  can  be  used  to  tune  the 
stability  by  using  divalent  cations  of  different  sizes.  Here  we  consid¬ 
ered  Sr-*Ca  substitution  and  A-cation  mixing,  to  give  the  following 
compounds:  LaSr2Cu2Ga07  -*■  LaCa2Cu2Ga07.  Without  A-cation 
mixing,  as  shown  in  Figure  10a,  we  found  a  strong  competition  be¬ 
tween  the  two  phases  with  ~  3  meV  favoring  the  NCSM.  However,  dif¬ 
ferently  from  LSCGO,  this  configuration  is  energetically  disfavoured 
(-680  meV)  with  respect  to  the  situation  with  A-cation  mixing  (Fig¬ 
ure  10b).  Importantly,  A-cation  mixing  stabilizes  the  Ima2  phase 
of  LaCa2Cu2Ga07  over  the  Pnma  one  by  -25  meV.  The  difference 
in  stability  between  LaSr2Cu2Ga07  and  LaCa2Cu2Ga07  can  be  ex¬ 
plained  considering  the  size  of  the  cations  involved  in  the  substitution. 

Ca2+  is  smaller  than  Sr2+  and  prefers  A-sites  in  layered  cuprates  such 


679.5  - 


679 


g  678.5  - 


678  - 


677.5 


1 

30 

1 

- 

25 

- 

- 

- 

20 

- 

- 

- 

15 

- 

- 

- 

10 

5 

0 

i 

- 

Ima2  Pnma 


Ima2  Pnma 

Figure  10:  Energy  difference  between  the 
Ima2  and  Pnma  structure  for  LaCajCuiGaCb 
without  (a)  and  with  (b)  A-cation  mixing. 


as  the  8-coordinate  site  between  cuprate  planes,  as  in  HgBa2CaCu206+a.~  The  chemical  substitution  reduces 
the  interlayer  space  and  demonstrates  that  chemical  pressure  is  a  promising  way  to  achieve  a  polar  conducting 
cuprate. 

Next,  based  on  under¬ 
standing  the  role  of  the  in¬ 
terlayer  Ga04  interactions 
(Figure  9),  we  identified  an¬ 
other  convenient  approach 
to  reduce  the  real-space  dis¬ 
tance  between  sub-blocks. 

Here  we  leveraged  the  out- 
of-plane  lattice  response  of 
a  thin  film  under  tensile 
epitaxial  strain  to  tune  the 
phase  stability.  Six  different 
substrates  were  surveyed 
in  the  range  bounded  by 

LaAlC>3  (3.821  A)  and  PrScC>3  (4.013  A).  Figure  11  shows  the  results  of  our  DFT  calculations  on  the  ener¬ 
getic  stability  of  the  polar  (Ima2)  phase  relative  to  the  centrosymmetric  (Pnma)  structure  for  thin  films 
of  LnSr2Cu2Ga07  (Ln=La  and  Y).  Note  that  in  the  compound  Yi_xCaxSr2GaCu207±5  superconductivity 
was  observed^  for  the  range  of  Ca-content  x  >  0.4.  We  therefore  also  computed  the  stability  of  ordered 
Y1/2Ca1/2Sr2GaCu207  compound.  We  find  that  epitaxial  strain  lifts  the  energetic  degeneracy  between  the 
phases.  Consistent  with  our  previous  discussion,  increasing  tensile  strain  provides  additional  stability  for 
the  polar  structure  through  controlled  variation  in  the  interlayer  d  spacing.  We  therefore  predict  that  these 
layered  Ga-containing  cuprates  when  grown  under  tensile  strain  should  be  polar  metals. 

Over  the  program  we  actively  engaged  and  continue  to  engage  two  experimental  groups  to  demonstrate  the 
phase  stability  of  the  polar  cuprate  in  bulk  and  thin  film  form.  First,  we  collaborated  with  Prof.  Ken  Poeppelmeier 
in  the  Chemistry  Department  at  Northwestern  University.  His  group  synthesized  ceramic  targets  for  pulsed- 
laser  deposition  (PLD)  growth  of  LaSr2Cu2Ga07  and  Yi_xCaA:Sr2Cu2Ga07  with  x  =  0.0, 0.2, 0.3, 0.4, 0.5.  For 
0  <  x  <  0.3  the  ceramics  can  be  made  phase  pure  with  small  particulate  size.  However,  for  x  >  0.4  it  is  much 
more  difficult  to  incorporate  Ca  into  the  lattice  (Figure  12). 


Figure  1 1 :  Effect  of  epitaxial  strain  on  the  energetic  stability  of  the  polar  Imal  phase  in 
LnSr2Cu2Ga07.  The  polar  structure  is  more  stable  (negative)  with  increasing  epitaxial  strain 
and  for  smaller  lanthanide  (Ln=Y)  cations.  Potential  substrate  materials  are  indicated.  The 
behavior  of  dc„  is  consistent  with  our  model  for  the  Ga04  interactions. 
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Using  the  prediction  of  large  tensile  strain  as  a 
route  to  stablize  the  polar  cuprate  phases,  the  Poep- 
pelmeier  group  also  grew  single  crystals  of  PrSc03 
as  a  substrate  material,  for  which  there  are  no  com¬ 
mercial  suppliers,  for  the  PLD  growth  by  Prof.  Yuri 
Suzuki’s  group  in  the  Applied  Physics  Department 
at  Stanford  University.  The  scandate  single  crystals 
were  grown  in  an  optical  float  zone  furnace  from  which  substrates  were  cut  and  Laue  diffraction  confirms  the 
zone  axes  of  the  single  crystals.  At  present,  the  Suzuki  group  is  determining  the  optimal  pulse-laser  deposition 
(PLD)  growth  conditions  on  LAO,  STO,  and  LSAT  using  the  ceramic  targets  prepared  by  the  Northwestern 
group.  We  are  hopeful  that  thin  films  with  the  correct  stoichiometry  will  be  realized. 


Figure  1 2:  Sintered  ceramic  and  PrSc03  single  crystals  by  floating- 
zone  method  and  corresponding  scandate  Laue  diffraction  pattern. 


B.4.  Polar  Metals  as  a  Platform  for  Multiferroic  Discovery. 

The  main  objective  of  this  effort  was  to  design  a  Mott-induced  multiferroic  by  taking  a  non-magnetic  polar 
metal  close  to  a  Mott  transition  and  inducing  a  magnetic  insulating  state  through  heterostructuring.  The 
motivation  was  to  realize  a  new  type  of  magnetoelectric  material,  i.  e. ,  where  their  is  strong  coupling  between  the 
primary  electric  and  magnetic  polarizations  owing  to  the  coexistence  of  ferroelectricity  and  ferromagnetism, 
from  a  NCSM.  Such  materials  are  highly  sought  for  in  low-power  multi-state  ferroelectric  random  access 
memories  (FE-RAM). 

We  examined  the  layered  ordered  (Li0s03)i/(LiNb03)i  super¬ 
lattice  formed  by  interleaving  LiNb03  and  LiOsC>3,  which  both  have 
polar  crystal  structures  (R3c).  '  The  former  is  an  insulator  with 
a  band  gap  of  3.25  eV  (Figure  13a)  whereas  Li0s03  is  metallic  with 
the  Fermi  level  dominated  by  the  Os  5 d- states  (Figure  13b).  First,  we 
examined  the  electronic  properties  of  the  newly  discovered  “ferro¬ 
electric  metal”  Li0s03  combining  density-functional  and  dynamical 
mean-field  theories.  We  showed  that  the  material  is  close  to  a  Mott 
transition.  We  identified  for  the  first  time  that  the  Hund’s  cou¬ 
pling  in  LiOsOs  make  it  an  ideal  candidate  for  realizing  a  Mott 
MF  due  to  the  multi-orbital  t2g  physics. 

Next,  we  describe  the  design  of  a  new  multiferroic  by  control 
of  the  electronic  structure  through  atomic  scale  engineering  of 
a  Mott  metal-insulator  transition  (MIT)  in  an  ultrashort  period 
(Li0s03)i/(LiNb03)i  superlattice.  We  used  electronic  structure  cal¬ 
culations  to  predict  that  the  (Li0s03)i/(LiNb03)i  superlattice  ex¬ 
hibits  strong  coupling  between  magnetic  and  ferroelectric  degrees 
of  freedom  with  a  ferrolectric  polarization  of  41.2  uC  cm~2,  Curie 
temperature  of  927  K,  and  Neel  temperature  of  379  K.  The  ferroelec- 


Figure  1 3:  (a)  The  Li0s03/LiNb03  superlattice 
exhibits  the  a~b~b~  tilt  pattern.  Atom-  and 
orbital-resolved  DOS  for  (b)  LiOs03,  (c)  LiNb03 
and  (d)  Li0s03/LiNb03  at  the  DFT-LDA  level. 


trie  properties  mainly  originate  from  cooperative  Li  and  O  displacements.  We  find  that  the  insulating  and 
magnetic  state  is  driven  by  an  enhancement  of  the  electronic  correlations  in  Li0s03  layers  owing  to  the  kinetic 
energy  reduction  of  the  t2g  orbitals  from  the  superlattice  geometry.  Our  efforts  have  uncovered  a  promising 
alternative  route  to  discovery  of  room-temperature  multiferroics:  One  could  search  for  correlated  polar 
metals  near  Mott  transitions  and  drive  the  phases  into  insulating  states,  rather  than  the  often-pursed  approach 
of  inducing  polar  displacements  in  robustly  insulating  magnets. 


B.5.  Correlation  and  Effects  in  Noncentrosymmetric  Hund’s  Metals. 

Using  a  first-principles  approach  based  on  density  functional  theory  and  dynamical  mean  field  theory,  we 
examined  the  electronic  properties  of  a  new  candidate  polar  metal  SrEiuVk^Of,  (SEMO).  Its  electronic  struc- 
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ture  shares  similarities  with  centrosymmetric  SrMo03  (SMO)  and  E11M0O3  (EMO),  from  which  it  may  be 
considered  an  ordered  derivative,  but  ferroelectric-like  distortions  of  the  divalent  cations  and  oxygen  anions 
lift  inversion  symmetry  mediated  by  an  anharmonic  lattice  interaction  in  the  metallic  state.  We  find  that 
Hund’s  coupling  promotes  the  effects  of  electronic  correlations  owing  to  the  Mo4+  d2  electronic  configuration, 
producing  a  correlated  metallic  phase  far  from  the  Mott  state  (Figure  14). 

The  ferroelectric-like  state  coexists  with  a  corre¬ 
lated  metallic  phase  which  relies  on  the  Elund’s  cou¬ 
pling  and  it  is  robust  with  respect  to  Mott  localization. 

The  contraindication  between  metallicity  and  polar  dis¬ 
tortions  is  thereby  alleviated  through  the  renormalized 
quasiparticles,  which  are  unable  to  fully  screen  the  or¬ 
dered  local  dipoles,  in  the  Hund’s  metal  SrEuMo20(v 

Our  main  finding  is  that  the  poorly  coherent  quasipar¬ 
ticles  are  ineffective  at  completely  screening  the  ordered 
local  dipole  moments.  These  results  strengthen  the  link 
between  correlated  metallic  states  and  the  propensity  of 
a  metallic  material  to  adopt  polar  structure,  facilitating 
the  selection  or  design  of  new  correlated  polar  met¬ 
als  for  enhanced  magnetoelectronic  responses  or  cus¬ 
tomized  antisymmetric  exchange  interactions,  which 
support  exotic  magnetic  textures  (helical  or  skyrmionic 
structures).  This  work  was  submitted  for  publication  to  Physical  Review  Letters,  and  can  be  found  online  as  a 
preprint  on  the  ArXiv. 

Although  our  current  understanding  of  polar  metals  reveals  that  most  of  them  are  in  a  correlated  elec¬ 
tronic  regime  with  bad  metallic  behavior,  we  do  not  eliminate  the  possibility  that  non-correlated  polar  metals 
may  exist;  rather,  strong  electronic  correlations  appears  to  be  a  favorable  ingredient  to  stabilizing  non- 
centrosymmetric  metallic  phases,  but  it  may  not  be  a  prerequisite  in  scenarios  where  the  inversion  lifting 
displacements  are  indeed  decoupled  from  the  low-energy  electronic  structure. 

B.6.  Polar  Distortions  in  Vanadates. 

Beyond  the  ruthenates  and  cuprates,  we  have  actively  been  exploring  the  structure-property  relationships  in 
other  transition  metal  oxides  which  are  close  the  the  metal-insulator  phase  boundary.  In  collaboration  with  the 
thin  film  growth  group  of  Prof.  Engel-Herbert  (Department  of  Materials  Science  and  Engineering,  PennState), 
we  investigated  the  strain  effect  on  the  electronic  and  magnetic  properties  of  the  correlated  metal  SrVC>3. 
The  Engel-Herbert  group  was  able  to  grow  oxygen  stoichiometric  SrV03  under  strain  which  exhibits  weak 
ferromagnetism.  We  note  that  in  the  literature  the  magnetic  behavior  of  SrVC>3  is  controversial.  Generally 
SrV03  is  described  as  a  Pauli  paramagnet.  However,  Dougier  et  al.  report  a  ferromagnetic  spin  density  wave 
in  bulk  SrVC>3  with  the  appearance  of  a  Pauli  paramagnetic  behavior  on  SrVC^-fl  for  S  >  0.1. 

In  Figure  15  we  show  the  calculated  magnetic  be-  20 

havior  of  SrVC>3  as  function  of  strain  within  Perdew-  > 

<D 

Burke-Ernzerhof  revised  for  solids  approximation  E_  10 

(PBEsol)  .  We  found  that  the  magnetic  behavior  of 
stoichiometric  SrV03  can  be  tuned  by  subtle  changes  0 

in  the  epitaxial  strain,  which  we  hypothesized  is  re¬ 
lated  to  the  Jahn-Teller  tendency  of  the  V4+  cation.  Figure  15:  Energetic  stability  of  various  magnetic  configurations 
The  ferromagnetic  phase  becomes  stable  for  small  of  SrV03  as  a  function  of  strain  calculated  with  PBEsol  functional, 
values  of  compressive  and  tensile  strain  consistent  with  the  experimental  observations.  Because  SrV03  is 
a  strongly  correlated  metal,  our  results  may  be  sensitive  to  correlations  effects  which  we  have  not  explicitly 


Figure  14:  (Color  online)  Quasiparticle  weights  (Z)  calculated 
within  DFT+DMFT  for  SMO,  EMO  and  SEMO  at  different  values 
of  the  ratio  Jh/U.  Inset:  Imaginary  part  of  the  inverse  dielectric 
function,  1/e,  for  SEMO  at  different  JH/U  values. 
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included.  For  this  reason,  we  performed  a  series  of  calculations  that  include  a  small  onsite  Coulomb  interac¬ 
tions  (17=1  eV).  We  found  essentially  the  same  results  for  small  values  of  strain.  Moreover,  we  analyzed  the 
possibility  to  induce  rotations  in  SrV03  as  function  of  strain  and  found  that  the  rotations  are  disfavored  in 
SrV03-  Note  that  as  part  of  this  collaboration  with  the  Engel-Hebert  group,  we  have  identified  that  ultrathin 
film  SrVOi  may  be  an  excellent  transparent-conducting  oxide. 

Finally,  we  examined  1/1  SrV03/CaV03  superlattices  (with  and  without  the  additional  Coulomb  inter¬ 
actions)  to  guide  our  search  for  probable  ground  state  atomic  configurations.  Our  main  results  are  that  (1) 
neglecting  the  correlations  produces  a  polar  and  metallic  ground  state  crystal  structure,  and  (2)  when  taking 
in  account  enhanced  correlations,  the  crystal  structure  becomes  centrosymmetric  with  space  group  Pi  and 
insulating  (U  >2  eV).  Such  a  correlation-induced  metal-insulator  transition  makes  the  superlattice  system, 
because  it  maybe  a  candidate  material  for  thermoelectric  applications. 

III.  Research  Training  and  Development 

The  program  supported  researchers  with  multiple  levels  of  expertise,  including  undergraduate  (A.  Cordi, 
who  received  a  M.S.  Degree)  and  graduate  (J.  Young)  students  and  post-doctoral  researchers  (D.  Puggioni 
and  M.  Gu).  These  researchers  have  become  proficient  or  fluent  in  performing  density  functional  theory 
calculations  with  the  Vienna  Ab  Initio  Simulation  Package  (VASP)  and  the  Ablnit  package  .  They  have  also 
learned  how  to  use  complimentary  post-processing  tools  to  understand  lattice  dynamical  simulations  and  aid 
in  the  application  of  group  theoretical  studies  to  explain  the  behavior  of  noncentrosymmetric  metals. 

To  provide  professional  development  opportunities,  all  researchers  supported  by  this  award  helped  to 
prepare  computational  supercomputing  hour  proposals,  from  which  the  results  described  before  were  obtained. 
These  efforts  include  the  submission  and  preparation  of  grant  proposals  to  the  Center  for  Nanoscale  Materials 
(CNM)  at  Argonne  National  Laboratory  and  the  Extreme  Science  and  Engineering  Discovery  Environment 
(XSEDE),  which  were  subsequently  awarded.  These  computational  hours  complement  those  made  available 
from  the  DOD  through  the  HPCMP.  All  researchers  also  presented  work  at  the  American  Physical  Society 
March  Meetings  in  2012-2015  and  at  other  conferences.  These  experiences  have  improved  each  researcher’s 
communication  skills  and  have  grown  each  persons  professional  network. 

IV.  Research  Dissemination 

Publications  The  award  supported  14  journal  article  publications.  See  the  ARO  Extranet  Report  for  details. 

Submitted  Manuscripts 

©  Lifting  Dirac  nodes  in  topological  semimetallic  perovskite  Srlr03  through  epitaxial  constraint,  J.  Liu, 
D.  Kriegner,  L.  Horak,  D.  Puggioni,  C.  Rayan  Serrao,  R.  Chen,  D.  Yi,  C.  Frontera,  V.  Holy,  A. 
Vishwanath,  J.M.  Rondinelli,  X.  Marti,  and  R.  Ramesh  (Submitted,  June  2015). 

©  Electronic  correlations  and  screening  effects  in  the  Hund’s polar  metal  SrEuMojO^,  G.  Giovannetti, 
D.  Puggioni,  J.M.  Rondinelli,  and  M.  Capone  (Submitted,  July  2015). 

Manuscripts  in  Preparation 

©  Origin  of  the  acentricity  of  the  layered  copper  oxide  superconductor  L  a  S  rj  Cui_  G  aOj,  D.  Puggioni, 
and  J.M.  Rondinelli 

©  Stabilization  of  the  polar  phases  in  LnSr2Cu2Ga07  (Ln-La,  Y),  D.  Puggioni,  and  J.M.  Rondinelli 
©  Polar  Metals:  Progress  and  Perspective,  Invited  Topical  Review  for  J.  Phys.  Condens.  Matter,  N.  A. 
Benedek,  and  J.M.  Rondinellit 

Presentations  See  electronic  entries  in  the  ARO  Extranet  Report. 
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